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ABSTRACT 

The ketenyl radical (HCCO) has recently been discovered in two cold dense clouds 
with a non-negligible abundance of a few 10 ” 11 (compared to H 2 ) (Agundez et al. 
2015). Until now, no chemical network has been able to reproduce this observation. 
We propose here a chemical scheme that can reproduce HCCO abundances together 
with HCO, H 2 CCO and CH 3 CHO in the dark clouds Lupus-lA and L486. The main 
formation pathway for HCCO is the OH + CCH —» HCCO + H reaction as suggested 
by Agundez et al. (2015) but with a much larger rate coefficient than used in current 
models. Since this reaction has never been studied experimentally or theoretically, this 
larger value is based on a comparison with other similar systems. 
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1 INTRODUCTION 

Cold dark clouds are assumed to be the most simple type of 
interstellar sources that can be used to test chemical mod¬ 
els before being applied to more complex sources such as 
protostars and protoplanetary disks (Bergin & Tafalla 2007; 
Agundez & Wakelam 2013). The clouds have characteristi¬ 
cally low temperatures (below 10 K), they are dense (be¬ 
tween a few 10 4 to a few 10 6 cm -3 ) and are shielded from 
any source of UV field. There exists however a source of 
UV photons produced by the emission from excited H 2 pro¬ 
duced by electrons originating from the ioniozation of H 2 by 
cosmic-ray particles (Prasad & Tarafdar 1983). More than 
70 molecular species (neutrals, cations and anions) have now 
been observed in these objects (Agundez & Wakelam 2013). 
HCCO was recently detected by Agundez et al. (2015) with 
a non negligible abundance of a few 10 -11 (compared to 
H 2 ); an abundance that they could not reproduce using their 
model. 

To study the formation and destruction of these species, 
chemical models have been developed considering an increas¬ 
ing number of processes both in the gas-phase and at the 
surface of the grains where a fraction of these species are 
believed to be formed (Herbst & van Dishoeck 2009). Such 
models are based on chemical networks containing chem¬ 
ical reactions with their associated rate coefficients. The 
networks have been constructed over time mostly driven 
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by new observations. Over the past few years, we have 
been revisiting those chemical networks from a chemists 
point of view, to provide a better description of the un¬ 
derlying processes. (Wakelam et al. 2009; Loison et al. 2012; 
Daranlot et al. 2012; Loison et al. 2014a,b). In this paper, 
we present a chemical scheme for the formation and destruc¬ 
tion of HCCO in cold dense clouds and compare our results 
with the observed abundances. 

2 CHEMICAL MODELING AND NETWORK 

2.1 The NAUTILUS chemical model 

The model used to simulate the abundance of HCCO in cold 
dense clouds is NAUTILUS. This model computes the abun¬ 
dance of molecules and atoms (neutrals but also some species 
in their cationic or anionic forms) in the gas-phase and also 
at the surface of interstellar grains. The equations and the 
chemical processes included in the model are described in 
earlier papers and we refer to Reboussin et al. (2014). 

For the gas-phase, many chemical reactions (including bi- 
nrolecular and unimolecular processes) are considered based 
on the kida.uva.2014 chemical network by Wakelam et al. 
(2015). Interstellar ice is modeled by a one-phase rate equa¬ 
tion approximation (Hasegawa et al. 1992): there is no dif¬ 
ferentiation between the species in the bulk and at the sur¬ 
face. The species from the gas-phase are allowed to physisorb 
on the surface of interstellar grains. Here, they can diffuse 
and react. The evaporation processes are: thermal (which are 
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inefficient at dense cloud temperatures), induced by cosmic- 
rays (following Hasegawa & Herbst 1993), and chemical 
desorption (as defined by Garrod et al. 2007). Photodes¬ 
orption has not been introduced in this model because 
photodesorption induced by cosmic-ray secondary photons 
was found to be inefficient compared to the chemical des¬ 
orption mechanism by Wakelam et al. (2014). In addition, 
Bertin et al. (2012) showed experimentally that photodes¬ 
orption of few CO layers (i.e. ~10) adsorbed on amorphous 
H 2 O was inefficient. Until more experiments are performed 
on water ices, it is reasonable to assume that molecular pho¬ 
todesorption is only a minor desorption process in dense 
clouds. Species can diffuse by thermal hopping only with 
a barrier of 0.5xEd (with Ed the species binding energy). 
Reactions with atomic hydrogen can undergo tunneling us¬ 
ing the formalism by Hasegawa et al. (1992). For chemical 
desorption, i.e. partial desorption of the products due to 
the exothermicity of the reactions occurring at the surface 
of the grains, approximately 1% of the products is allowed 
to desorb. In addition to the usual surface diffusion reac¬ 
tions, we have introduced the complexation and low temper¬ 
ature Eley-Rideal mechanism from Ruaud et al. (2015). The 
Cosmic-Ray Induced Diffusion mechanism is not included 
because it does not have any impact on the surface chem¬ 
istry at high visual extinction (see Reboussin et al. 2014). 
To simulate the chemistry of dense clouds, the model is used 
with homogeneous conditions and integrated over 10 7 yrs. 
The gas is initially composed of atoms (partly ionized de¬ 
pending on the ionization potential of the atoms). The abun¬ 
dances are the same as in Hincelin et al. (2011). For fluorine, 
not included in Hincelin et al., we use the depleted value of 
6.68 x 10 -9 (compared to the total proton density) from 
Neufeld et al. (2005). The C/O elemental ratio is the sub¬ 
ject of considerable debate and is of crucial importance for 
the modeling of dense clouds (Hincelin et al. 2011). For our 
standard model, we have used a C/O ratio of 0.7 (i.e. the 
oxygen elemental abundance is 2.4 x 10~ 4 ). The model was 
run with a dust and gas temperature of 10 K, a total pro¬ 
ton density of 2 x 10 4 cm -3 , a cosmic-ray ionization rate of 
1.3 x 10 -17 s _1 , and a visual extinction of 30. 


2.2 Modification of the network 

The initial gas-phase network is kida.uva.2014 
Wakelam et al. (2015). The network for surface reac¬ 
tions and gas-grain interactions is based on the one from 
Garrod et al. (2007) with several additional processes from 
Ruaud et al. (2015). The HCCO molecule was introduced 
by Ruaud et al. but was not discussed in their paper 
since it was mostly an intermediate to the formation of 
more complex species on the surface. Indeed, the gas-phase 
reactions involving this species were very limited and as we 
will show in the next section, the gas-phase abundance of 
this species was very low at times larger than 10 s yr. In 
addition to these reactions (listed in Tables A.l and A.2 
of Ruaud et al. 2015), we have introduced and reviewed 
reactions involving HCCO radicals (presented in Table Al) 
to take into account the most important pathways for the 
formation and destruction of this species in the gas-phase. 



Figure 1. HCCO abundance (with respect to H 2 ) predicted by 
our model for typical dense cloud conditions (see section 2.1) as a 
function of time in the gas-phase (solid lines) and at the surface 
of the grains (dashed lines). Grey lines show the results obtained 
with the previous network while the black lines show the results 
including the new chemistry. The light grey curve for HCCO in 
the gas phase has been obtained with the new chemistry but 
removing the reaction s-CCO + s-H on the surface. 

3 PREDICTED ABUNDANCES 

Using the model described above, we have run the simu¬ 
lations with the new chemistry and compared the predic¬ 
tions with the previous results. The HCCO abundance pre¬ 
dicted by the model in the gas-phase and at the surface 
of the grains is shown in Fig. 1. The formation of HCCO 
on the surface is very efficient between 10 2 and 10 s yr be¬ 
cause of the reaction s — H + s — CCO —> s — HCCO, 
which is barrier-less for the equivalent gas phase process 
(Bauer et al. 1985; Schmidt et al. 1983; Horie et al. 1983). 
s-CCO is formed through the Eley-Rideal mechanism devel¬ 
oped in Ruaud et al. (2015): gas-phase carbon atoms react 
without a barrier with CO on the ice surface. The absence 
of a barrier for the C + s-CO —» s-CCO reaction is deduced 
from the study of Husain & Kirsch (1971) and ab-initio cal¬ 
culations performed by Ruaud et al. (2015). The gas-phase 
abundance of HCCO before 10 5 yr is then the result of the 
chemical desorption (see section 2.1) of s-HCCO during its 
formation on the surface. The abundance of HCCO in the 
gas-phase predicted by the new model without this reaction 
is also shown in Fig. 1. Note that in Table B1 of Ruaud et al. 
(2015) only the channel leading to s-HCCO is indicated but 
the chemical desorption channel was also included. 

After 10 5 yr, its abundance decreases in the previous 
model because s-HCCO is also hydrogenated to give s- 
H 2 CCO and we did not have any efficient gas-phase 
formation processes. With the new chemistry, its abun¬ 
dance in the gas-phase increases at larger times. The 
main production reaction is OH + CCH —>- HCCO + H 
as already suggested by Agundez et al. (2015). There are 
no theoretical or experimental studies of this reaction. 
Agundez et al. used a rate coefficient of 3 x 10 _n cm 3 s _1 
based on Frenklach et al. (1992). This value is however 
simply an estimation and there is no justification in 
Frenklach et al. (1992) of this choice. Based on similar 
systems, i.e. the O + CCH rate constant (Boullart et al. 
1996; Devriendt et al. 1996; Devriendt & Peteers 1997; 
Georgievskii & Klippenstein 2011) or OH + radicals re¬ 
actions (De Avillez Pereira et al. 1997; Jasper et al. 2007; 
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Sangwan et al. 2010), we propose the use of a larger rate 
coefficient of 2 x 10 ~ 10 cm 3 s -1 . This rate constant is close to 
the capture rate dominated at low temperature by dipole- 
dipole interactions and to the one precisely computed by 
Georgievskii & Klippenstein (2011) for O + CCH. The rate 
constant calculation methodology has been described in de¬ 
tail in Loison et al. (2014a). 

In this case, the formation of HCCO is faster than 
its destruction through reactions with H, O and Hg\ 
There is also another potential source of HCCO radicals: 
the CCH + O 2 reaction. The rate constant of this re¬ 
action is high at low temperature (Chastaing et al. 1998; 
Vakhtin et al. 2001) and HCCO has been experimentally 
identified as one of the products (Lange & Gg. Wagner 
1975) with a non-negligible branching ratio between 20% 
to 40% (Lange & Gg. Wagner 1975; Laufer & Lechleider 
1984). The production of HCCO radicals is also suggested 
from theoretical calculations (Sumathi et al. 1998; Li et al. 
2004). However in the Lange & Gg. Wagner (1975) and 
Laufer & Lechleider (1984) studies, the rate constant of the 
CCH + O 2 reaction was underestimated and the rate con¬ 
stant of the HCCO + O 2 reaction was overestimated leading 
to a large uncertainty in the HCCO product branching ratio. 
Moreover our current astrochemical model leads to an O 2 
abundance higher than observed for dense molecular clouds 
(Pagani et al. 2003), inducing an overestimation of HCCO 
production through the CCH + O 2 reaction. Instead of in¬ 
cluding HCCO production from the CCH + O 2 reaction in 
our standard model (the CCH + O 2 reaction producing only 
H, HCO, CO and CO 2 Su et al. 2000), we performed a test 
run to see the effect of its inclusion. When the CCH + O 2 —» 
HCCO + O reaction is included, an increase of the HCCO 
abundance is observed in the 3-5 x 10 s yr range (reaching few 
10 -11 of H 2 ) corresponding to the peak of the O 2 abundance 
profile. 

Fig. 2 shows the abundances of HCO, H 2 CCO and CH 3 CHO 
predicted by our new model as a function of time for the 
physical conditions described in section 2.1. Compared to 
Agundez et al. (2015), H 2 CCO is predicted to be much less 
abundant very likely because we have introduced its de¬ 
struction reaction with atomic carbon with a rate coeffi¬ 
cient of 3 x 10 -1 ° cm 3 s _1 (see Ruaud et al. 2015). Our 
predicted CH 3 CHO abundance is quite different from the 
one of Agundez et al. (2015) with the UMIST RATE12 or 
the kida.uva.2014 networks. The difference with their re¬ 
sults using the kida.uva.2014 network is mainly due to the 
fact that they do not consider reactions at the surface of the 
grains. In our present simulation, CH 3 CHO is produced in 
the gas phase through the O + C 2 H 5 reaction, C 2 H 5 being 
efficiently produced through surface reactions (mainly s-H 
+ S-C 2 H 4 —»• C 2 H 5 ). The difference between the gas phase 
models of Agundez et al. (2015) using the two networks is 
more complex to analyse. One main difference though be¬ 
tween our network and UMIST RATE12 is the inclusion of 
destruction reactions with atomic carbon. 


4 COMPARISON WITH OBSERVATIONS 

HCCO has only been observed in two dark clouds, Lupus- 
1A and L483 by Agundez et al. (2015) with an abundance 
of a few 10 -11 compared to H 2 . These authors also give 



Figure 2. HCO, H 2 CCO and CH 3 CHO abundances (with re¬ 
spect to H 2 ) predicted by our model using the new chemistry for 
the typical dense cloud conditions described in section 2.1 as a 
function of time in the gas-phase. 


upper limits in three other cold sources with a very strong 
limit for L1521F of 7 x 10 -13 . They also detected HCO, 
H 2 CCO and CH 3 CHO in the same sources. We will compare 
our model results with the abundances listed in Table 2 of 
Agundez et al. (2015) keeping in mind the uncertainties in 
the H 2 column densities in these sources (see Agundez et al. 
2015, for discussion). In TMC-1, HCCO was not detected 
and an upper limit of 2.9 x 10 ~ 12 (compared to H 2 ) was de¬ 
rived by Turner & Sears (1989). The reanalysis of the HCCO 
observation by Agundez et al. (2015) gave an upper limit of 
5 x 10 -11 . 

For all the molecules (Figs. 1 and 2), we are able to reproduce 
reasonably well the observed abundances for times between 
10 5 and 10 ® yr (chemical ages usually assumed for dark 
clouds). Only HCO is slightly overproduced by our model. 
It is interesting to note that the predicted HCO abundance 
presents a smooth increase while H 2 CCO and in particular 
CH 3 CHO have peak abundances around 2 x 10 5 yr. This 
peak abundance of H 2 CCO is ten times larger than the ob¬ 
served one in any source. The strong increase of H 2 CCO and 
CH 3 CHO near 10 5 yr is due to the fact that these closed shell 
molecules react only with C atoms and not with H, O or N 
atoms; the carbon atom abundance strongly decreasing be¬ 
tween 10 5 and 2 x 10 5 yr. The relatively stable abundance of 
HCO with time may explain the uniform abundance of this 
molecule observed by Agundez et al. (2015) in all sources. 
The main production and destruction reactions for those 
molecules are summarized in Table Bl. 

The observed sources present some variations in density 
and temperature that may influence the chemistry (see 
Agundez et al. 2015, and references therein). To test this, we 
have run our model with a larger total proton density (nH = 
n(H) + n(H 2 )) of 2 x 10 5 and even 2 x 10 6 cm -3 (see Fig. 3), 
keeping the other parameters similar to the ones described 
in section 2 . 1 . As a general result of more collisions, the pre¬ 
dicted gas-phase abundances are increased earlier but the 
peak abundances and late time abundances are smaller be¬ 
cause of more efficient depletion on the grains. If some of the 
sources are denser such as Lupus-1A and L1521F, it is very 
likely that they have not been that dense for a long period of 
time. Some of the sources may also be slightly warmer than 
the typical 10 K usually assumed. We have run our model 
with gas and grain temperatures of 15 and 20 K (see Fig. 3). 
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Figure 3. HCO, HCCO, H 2 CCO and CH 3 CHO abundances 
(with respect to H 2 ) predicted by our model using the new chem¬ 
istry as a function of time in the gas-phase. Some of the phys¬ 
ical parameters have been changed compared to the standard 
model: the gas and dust temperature is set to 15 K for the solid 
lines and to 20 K for the dashed lines (keeping the H density to 
2 X 10 4 cm -3 ) ; the density was set to 2 X 10 s cm -3 for the dotted- 
dashed lines and to 2 X 10® cm' 3 for the dotted lines (keeping 
the temperature to 10 K). 


The 15 K model is not very different from the 10 K standard 
model. In contrast, the 20 K model produces fewer molecules 
in general except for CH 3 CHO which remains relatively un¬ 
affected. After 10 5 yr, the 20 K regime is characterized by a 
smaller abundance of CCH and C 2 H 3 , a smaller abundance 
of CO onto the grain surfaces (the evaporation temperature 
of CO being around 18 K), a larger abundance of atomic hy¬ 
drogen in the gas-phase, and an increase of the diffusion rate 
of heavy radicals onto the surface. Before 10 s yr, less HCCO 
is produced because of a smaller abundance of CO on the 
grains (less CCO is formed). After 10 5 yr, it is the larger gas- 
phase abundance of atomic hydrogen that destroys HCCO 
efficiently in the gas-phase. The lower abundance of HCCO 
in the ices leads to less H 2 CCO in the gas-phase. 

Using an approach similar to Loison et al. (2014a) to com¬ 
pare the HCO, HCCO, H 2 CCO and CH 3 CHO observed 
and modeled abundances obtained with the standard phys¬ 
ical model, we have found that the observed abundances 
in Lupus-IA and L483 are best reproduced at 1.8 x 10 s yr 
with a mean difference between observed and modeled abun¬ 
dances of less than a factor of 3. Using a gas and dust tem¬ 
perature of 15 K, the best mean agreement is smaller than 
a factor of 2 at 1.1 x 10 s yr because HCCO and H 2 CCO 
are better reproduced. Considering both the uncertainties 
in the observations and the chemical modeling, this differ¬ 
ence between the models is however insignificant. The ’’best 
age” that we give here is only based on the comparison be¬ 
tween the modeled and observed abundances of the four 
molecules. In addition, estimating an age for such sources 
requires us to define the time zero as a starting point. The 
formation of dense clouds is a continuous process from the 
diffuse medium. The formation of molecules is likely to start 
early on during this process when the density and the vi¬ 
sual extinction are smaller than the ones considered in our 
model. For these reasons, this age should only be considered 
as the best time in our model to reproduce the observations. 
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Table Al. List of gas-phase reactions added to the model and associated parameters. 



Reaction 


A E (kJ/mol) 

a 

P 

7 

Fo 

g 

Reference 

1 

H + HCCO 

-+ ch 2 + CO 

-113 

1.7 x 10~ 10 

0.17 

0 

1.6 

0 

1,2 

2 

C+ + HCCO 

-+ C 2 H+ + CO 

-411 

2.0 x 10 -9 

0 

0 

3 

0 

3 



-+ HCCO+ + C 

-143 

2.0 x 10" 9 

0 

0 

3 

0 

3 

3 

C + HCCO 

-+ CCH + CO 

-437 

2.0 x 10~ 10 

0 

0 

3 

0 

4 

4 

N + HCCO 

-+ HCN + CO 

-624 

2.0 x 10 -11 

0.17 

0 

3 

10 

5 



HNC + CO 

-569 

1.0 x 10“ u 

0.17 

0 

3 

10 

5 



> NCCO 1 H 

-220 

1.0 x KT 11 

©++ 

0 

■3 

4& 

5 

5 

O + HCCO 

-+ H + CO + CO 

-428 

1.6 x 10“ 10 

0 

0 

1.6 

0 

2 



-+ CH + C0 2 

-221 

4.9 x 10" 11 

0 

560 

2 

100 

2 

6 

OH + CCH 

-+ HCCO + H 

-220 

2.0 x 10“ 10 

0 

0 

3 

0 

3 



-+ CH + CO + H 

+96 

0 

0 

0 

3 

0 

3 

7 

HCCO + H+ 

-+ H 2 CCO+ + H 2 

-387 

1.0 

3.1 x 10" 9 

2.8 

2 

0 

6 

8 

HCCO + H 3 0+ 

-+ H 2 CCO+ + h 2 o 

-122 

1.0 

1.4 x 10“ 9 

2.8 

2 

0 

6 

9 

HCCO + HCO+ 

-+ h 2 cco+ + CO 

-244 

1.0 

1.3 x 10' 9 

2.8 

2 

0 

6 

10 

HCCO + N 2 H+ 

-+ h 2 cco+ + n 2 

-339 

1.0 

1.3 x 10~ 9 

2.8 

2 

0 

6 

11 

H 2 CCO+ + e" 

-+ ch 2 + CO 

-607 

2.0 x 10“ 7 

-0.5 

0 

3 

0 




—> H + HCCO 

-495 

1.0 x 10“ 8 

-0.5 

0 

3 

0 




-> H + CH + CO 

-178 

1.0 x 10“ 7 

-0.5 

0 

3 

0 




—s> H + H + CCO 

-92 

0 

0 

0 

3 

0 


12 

CH3CO+ + e" 

-+ ch 3 + CO 

-629 

1.0 x 10“ 7 

-0.5 

0 

3 

0 




-> H + H 2 CCO 

-509 

1.0 x 10“ 7 

-0.5 

0 

3 

0 




—s> H + CH 2 + CO 

-185 

1.0 x 10“ 7 

-0.5 

0 

3 

0 




-+ H + H + HCCO 

-72 

0 

0 

0 

3 

0 



A E is the enthalpy of reaction. 

a, (3 ,7 are the parameters to compute the temperature dependent rate coefficients. The Arrhenius-Kooij formula described in 
Wakelam et al. (2012) should be used for the reactions except for reactions 7 to 10 for which the ionpoll formula from the KIDA 
database (http://kida.obs.u-bordeauxl.fr) should be used (see also Wakelam et al. 2012). The rates can be used for temperatures 

between 10 and 300 K. 

Fo and g are parameters used to describe the temperature dependent uncertainty in the rate coefficient (see Wakelam et al. 2012). 

References: 

1 Glass et al. (2000) 

2 Baulch et al. (2005) 

3 Close to capture rate 

4 Estimation considering the well known reactivity of carbon atoms 
5 Reaction on triplet surface with some spin orbit crossing. Part of the HCN should lead to HNC considering the available energy. The 

last channel was removed because NCCO is not in KIDA. 

6 Capture rate 
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Table Bl. Main reactions of production and destruction for HCO, HCCO, H 2 CCO and CH 3 CHO for two different times and using the 
new model. 


Species 

10 4 yr 

Formation 

Destruction 

HCO 

O + CH 2 -4 H + HCO 

O + C 2 H+ -4 HCO + H+ 

H + HCO -4 CO + H2 

C 4- HCO -4 CH + CO / H + CCO 

HCCO 

s-H + s-CCO -4 HCCO 

O + HCCO -4 H + CO + CO 

C + HCCO -4 CCH + CO 

H2CCO 

O + C2H3 -4 H2CCO + H 
s-H + s-HCCO -4 H 2 CCO 

C + H2CCO -4 C2H2 + CO 

CH3CHO 

O + C2H5 —> H + CH3CHO 

C + CH3CHO -4 C2H4 + CO 

Species 

10 6 yr 

Formation 

Destruction 

HCO 

s-H + s-CO -4 HCO 
H2COH+ + e“ -4 H + H + HCO 

H + HCO -4 CO + H 2 

HCCO 

OH + CCH -4 HCCO + H 

H + HCCO -4 CH 2 + CO 

H2CCO 

O + C2H3 -4 H + H2CCO 

H+ + H2CCO -4 H 2 + CH3CHOH+ 

CH3CHO 

s-H + S-CH3CO -4 CH3CHO 

H+ + CH3CHO -4 H 2 + CH3CHOH+ 
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